Existence and merging of Dirac points in a-(BEDT-TTF) 2 I 3 organic conductor 
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We reexamine the existence and stability conditions of Dirac points between valence and conduc- 
tion bands of 3/4 filled a-(BEDT-TTF)2l3 conducting plane. We consider the usual nearest neigbhor 
tight binding model with the seven transfer energies that depend on the applied pressure. Owing 
to the four distinct molecules (A, A', B, C) per unit cell of the Bravais lattice, the corresponding 
Bloch Hamiltonian is a 4 x 4 matrix T-i(k) for each wave vector k of the Brillouin zone. In most 
previous works the study of Dirac points was achieved through direct numerical diagonalization 
of matrix T-i{k). In this work we develop a novel analytical approach which allows to analyze the 
existence and stability conditions of Dirac points from the knowledge of their merging properties 
at time reversal points. Within this approach we can discuss thoroughly the role of each transfer 
integrals and demonstrate that inversion symmetry is convenient but not necessary. 



I. INTRODUCTION 

The anomalous temperature dependence of Hall co- 
efficient in organic conductor a-(BEDT-TTF) 2 I 3 under 
pressure 1,2 has been an unresolved problem for a long 
time. The issue was recently clarified by tight-biding 
band structure calculations which point out that the ori- 
gin of this anomalous behaviour is the existence of a 
zero gap state (ZGS) due to the presence of Dirac points 
(between valence and conduction bands) exactly at the 
Fermi energy 3 . Indeed when Dirac points appear at the 
Fermi energy, then electrons and holes pockets eventu- 
ally disappear resulting in a semimctal phase with a van- 
ishing density of states at the Fermi energy. Further 
tight-binding calculations, based on diffraction data un- 
der pressure 4 , have established the robustness of the ZGS 
under pressure. This robustness relies on the stability of 
Dirac points at the Fermi energy which itself is attributed 
to the persistence of an inversion symmetry under pres- 
sure. 

The scenario of the ZGS has allowed much progress 
in the understanding of the physical properties of a- 
(BEDT-TTF) 2 I 3 5 ~ 7 . Nevertheless the conditions for the 
existence and stability of Dirac points at the Fermi en- 
ergy are still not well understood. Indeed compared to 
graphene that has only two atoms per unit cell, the sim- 
plest tight-binding model of a-(BEDT-TTF) 2 I 3 conduct- 
ing plane contains six electrons for four distinct molecules 
A, A' , B, C per unit cells and seven distinct nearest neigh- 
bor transfer energies between them. The Bloch Hamil- 
tonian matrix U(k) associated to a-(BEDT-TTF) 2 I 3 is 
then of size 4x4 thus preventing much analytical treat- 
ment as opposed to tight-binding model of graphene that 
leads to a 2 x 2 Bloch hamiltonian matrix which allows full 
analytical understanding of Dirac point properties 8 ' 9 . As 
a consequence, much understanding of the existence and 
stability of Dirac points in 3/4 filled a-(BEDT-TTF) 2 I 3 
conducting plane relies on numerical diagonalization of 
the 4x4 Bloch Hamiltonian matrix H(k). 

Many numerical studies have clearly established the 



existence and stability of Dirac points under varying pres- 
sure (e.g. pressure effect is encoded in the value taken 
by the seven transfer energies). They have shown that 
Dirac points move in the Brillouin zone by varying pres- 
sure and can merge at time reversal points for some crit- 
ical pressure such that the system is fully gapped above 
this critical pressure. Further studies have also exam- 
ined the role of onsite potentials either due to anions or 
Hartree mean-field electron-electron interaction. As a re- 
sult it was shown that when such potentials break inver- 
sion symmetry between molecules A and A' this prevents 
the appearance of Dirac points 10-12 

Apart from these numerics, recent analytical ap- 
proaches have obtained interesting results. On the one 
hand, by neglecting some transfer energies one can built 
an effective 2x2 Bloch Hamiltonian matrix with only 
four transfer energies. It has been shown that this ef- 
fective model leads to the so called tilted Dirac cones 
with very anisotropic velocities 3 . More recently an exact 
mapping of the 4x4 Hamiltonian matrix to a 2 x 2 effec- 
tive matrix model was built 13 . This 2x2 effective model 
encodes the Dirac points physics in a self-consistent man- 
ner, it is however not clear how to reconstruct from such a 
reduced 2x2 matrix the physical properties that concern 
only the valence and conduction bands of the full 4x4 
Hamiltonian matrix. On the other hand Mori 14 recently 
examined a simplified 4x4 Bloch Hamiltonian matrix, 
constrained by inversion symmetry, in which he consid- 
ers only interchains transfer energies t\, n (n = 1,2,3,4) 
and where transfer energies t an (n = 1,2,3) along the 
stacking axis are neglected such that no direct A — A' 
and B — C hopping are allowed. This model appears ex- 
actly solvable and allows full analytical understanding of 
the existence and stability conditions of Dirac points and 
their k space motion under pressure. 

Despite all these numerical and analytical studies, the 
question of the origin of the stability of the Dirac points 
in 3/4 filled a-(BEDT-TTF) 2 I 3 conducting plane is still 
pending. Is inversion symmetry necessary or only con- 
venient ? How important is the sign and modulus of 
each transfer energies for the stability of Dirac points ? 



Answering these questions remains crucial, because once 
the Dirac points are well established it is then possible to 
consider only the 2x2 effective low energy model that de- 
scribes the valence and conduction bands in the neighbor- 
hood of the Dirac point. Such low energy model allows 
then to examine the novel physical properties associated 
to the presence of ZGS with tilted Dirac cones 10,11,15,16 . 

The aim of the present paper is to present a novel 
method to examine the existence and stability of Dirac 
points in 3/4 filled a-(BEDT-TTF) 2 I 3 conducting plane. 
Shortly said, our method allows to obtain analytical al- 
beit complicated expressions of eigen energy bands of the 
full 4x4 Bloch Hamiltonian matrix H(k) in the presence 
of the seven transfer energies tbn and t an . However we 
show that to characterized the properties of Dirac points 
between valence and conduction bands it is sufficient to 
study the properties of an alternative quantity K(k) that 
has a simpler analytical expression. We stress that equiv- 
alent quantities could be defined to study contact points 
between other neighboring bands as it was exemplified in 
a recent numerical study 17 . Using this quantity K(k) 
we can then examine thoroughly the role played by the 
different transfer energies on the stability of Dirac points 
and discuss the importance of inversion symmetry. 

The paper is organized as follows. In section II we 
setup the notations and present in more details the tight- 
binding model and associated Bloch Hamiltonian matrix 
of 3/4 filled a-(BEDT-TTF) 2 I 3 conducting plane. We 
also present the different steps to dcrivate formal an- 
alytical expressions of eigen energy bands. In section 
III as a first application of the method, we reconsider 
the simplified case of Mori in which transfer energies t an 
{n = 1,2,3) along the stacking axis are neglected. We 
explain how to generalize Mori's result to system with- 
out inversion symmetry. We find that the absence of 
inversion symmetry is not detrimental to the existence 
and stability of Dirac points. A possible explanation of 
this unexpected stability is the existence of a chiral sym- 
metry which appears because the system becomes bipar- 
tite when the transfer energies t an are absent. In sec- 
tion IV we consider the full model with the seven trans- 
fer energies. We first introduce the alternative quantity 
K(k). Using this quantity K(k) we can then examine 
thoroughly the role played by the different transfer en- 
ergies on the existence and stability of Dirac points and 
discuss the importance of inversion symmetry. Section V 
gives a summary of our main results. 




FIG. 1. (Color online) Crystal structure on two-dimensional 
plane with four molecules A, A', B, and C, in the unit cell 
where the respective bonds represent the seven transfer ener- 
gies t a with a = al, • • ■ , £t>4- 



where a represents the four distinct molecules A, A', 
B, and C per cell and i(= Hi) represents the Bravais 
position of a given cell in the effective rectangular Bra- 
vais lattice. The quantity ^ denotes the electron cre- 
ation operator on molecule a in cell Rj of the Bravais 
lattice. Coefficients f a ,/3;i,j denote transfer energies be- 
tween nearest neighbor sites. As shown in Fig. I 18 there 
are seven transfer energies given by tbi ■ ■ ■ tbi along the 
direction of x-axis (b axis) and by t a i ■ ■ ■ i a3 along the 
direction of j/-axis ( stacking a-axis). A key properties 
of these a-organic material is that these transfer energies 
can vary with the applied pressure P. A linear interpo- 
lation based on diffraction data 4 under uniaxial pressure 



allows to write t a (P) = t°(l + K a P) = t° a (l+£-) where 



., ' " ^_ + 0n , p 

t Q a , K a and P a = 1/ K a are indicated in the following ta- 
ble I. Note that t a (-P a ) = and t a (P a ) = 2t° a . A small 
value of \P a \ indicates that the corresponding t a (P) is 
highly sensitive to pressure. The present choice of the 

a al a2 a3 bl b2 b3 b4 



t° [meV] 
K a [xHrtbar- 1 ; 
P a = K- 1 [kbar] 



-28 -48 20 123 140 62 25 
89 167 -25 11 32 
11.2 6 -40 - 91 31 - 



II. TIGHT-BINDING MODEL FOR 
a-(BEDT-TTF) 2 I 3 CONDUCTING PLANE 

The nearest neighbor tight binding Hamiltonian for a- 
(BEDT-TTF) 2 I 3 is given by 



i,j ct,/3 



(1) 



TABLE I. The seven transfer integrals t an ,tt, n depend linearly 
on the applied pressure P like t a (P) = t° a (1 + K a P) = t°(l + 
-p-) where t a and K a and the characteristic pressure P a are 
indicated for each bond. 



sign for each parameter t° a , K a is the same as that of 
the original one by Mori 18 , but is different from that of 
Katayama et al. 3 . 

For later use we also define t a ± = t a3 ± t a2 and sim- 
ilarly tb± — tb3 ± tbi and t c ± — tbi ± tbA with t c ± in- 



3 



dependent on pressure. Accordingly we rewrite each of 
them as t a (P) = t°(l + K a P) = t° a (l + p^) with the 
corresponding parameters t% , K a , P a written in table II. 
Anticipating on our results we point out three key prop- 



a+ 



b+ b- c+ c- 



t° a [meV] -28 68 202 -78 148 -98 

K a [xlO^kbar" 1 ] 304.1 110.5 17.4 -5.7 
P a = K' 1 [kbar] 3.3 9 57.5 -175.4 - - 

TABLE II. ta and K a and P a values of effective transfer en- 
ergies t a (P) = £°(1 + for a = (a±,b±,c±). 



\Un\ 



hrt' Pan 



< 1, (") 



> 1 



erties that can be read of table II: (i) 

^ p*" 1 '^ » 1 and more importantly (iii) ^ 

with f] = ±,t/ = ±. Lastly we also define the two ratio 

t°_ 

x± — |p— | that will appear many times in our calcu- 

b~ 

lations. Quantitatively with parameters of table II we 
obtain 



1.25 > 1, 



s + = |S±|=0.73<l. 



Bloch Hamiltonian matrix 



(2) 



We define Bloch state creation operators a' a (k) = 
V^ 1 J2j ciRjka aj, and four component operators a^(k) = 
(a^(fe), a\,(k), a^ B (k), a^(fe)). We can rewrite H = 
Ibz a){k)U{k)a{k) where the 4 x 4 matrix H(k) 
reads as 

( a b c\ 

a* d e 
b* d* / 
\c* e* f* 0/ 

with each matrix element given by 

a = t a3 + t a2 e lk v , 
b = t b3 + t b2 e lk * , 
c = tue ik " + t bl e lk * +ik * 



H(k) 



(3) 



d = t b2 + t b3 e tk * 
e = t bl + t bA & lk * 
f = t„i + t al e lk y 



(4) 



with —it < k x ,k y < 7T. A key property of the ma- 
trix elements is that the modulus of each of them has a 
monotonous dependency in either k x £ [0, 7r] for b, c, d, e 
or k y e [0, 7r] for a,/. As examples \b\ 2 



Xy 

l-x +2 
2 l b 

and 



l+x +2 



IS 



2 "6+ ^ 

cosk x (x e [1 : — 1]) 
a linear function of 



is a linear function of ir 

12 _ i+y+2 , l^J/ f 2 

y = cosfc y (y e [1 : -1]). 

B. Formal analytical derivation of eigen energy 
bands and gaps 

We denote E n (k), n — 1,2,3,4 the four eigen en- 
ergies (bands) obtained by diagonalizing matrix "H(fc) 
at a given k with E^k) > E 2 (k) > E 3 (k) > E 4 (k). 
Since a-(BEDT-TTF) 2 I 3 is 3/4 filled we call E 2 (k) va- 
lence band and Ei(k) conduction band even though there 
might be some band overlap. In this framework, by 
varying k over the whole Brillouin zone one can explore 
the possibility of Dirac point by searching for position 
fco around which the gap between valence band E 2 (k) 
and conduction band E\(k) vanishes linearly such that 
Ei(k a + q) — E 2 (k + q) = c q \q\. In most previous works, 
all this procedure is achieved through numerical diago- 
nalization of matrix H{k) and until now it has prevented 
a clear understanding of role of the different transfer en- 
ergies and of inversion symmetry on the existence and 
stability of Dirac points. 

We now explain how one can go further with analytics 
and derivate formal analytical albeit complicated expres- 
sions for energy bands. As we already pointed out in the 
introduction we are not really going to use these expres- 
sions, nevertheless the formal derivation appears quite 
instructive and moreover it allows us to reach equations 
(8) and (9) which are the starting points of the analysis 
of Dirac points properties presented in sections III and 
IV. 

To start with, we remind that E n {k) is a root of 
F k (uj) = T>et\H- uj\ with 



F fe (w) = lj a - B k uj 2 + C k u + D k , (5) 
where coefficients B k , C k ,D k are given by 



B k = \a\ 2 + \b 2 \ + \c\ 2 + \d\ 2 + \e\ 2 + \f\ 2 , 

Cfc = -(dfe* + ed*f* + a*e*c + c*ae + a*d*b + adb* + bfc* + b*f*c), 

Dk = |/| 2 M 2 + H 2 H 2 + |c| 2 |d| 2 - a*e*bf - aeb* f* - a*d*f*c- adfc* - b*e*cd - bec*d*. 



(6) 



Note that B k > for any k while the sign of C k and 
D k depends on k. The matrix Ji k is traceless and then 
we have J2 n E n (k) = 0. Owing to this property the first 



step consists to rewrite energy bands E n {k) as 

Et(k) = E(k) + A(k), 
E 2 (k) = E(k) - A(fe), 
E 3 (k) = -E(k)+A'(k), 
E 4 (k) = -E(k)-A'(k), 



(7) 
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with E, A, A' > and A + A' < 2E for any k. For a 
3/4 filled system, E(k) is the energy value located in the 
middle of the gap 2A(fc) separating valence and conduc- 
tion bands. The second step consists to relate E, A, A' 
to the three coefficients Bk, Ck, Dk'- 



(n = 1, 2, 3). In that situation a 
Bk 

c k = o, 



f = and 



\b 2 \ + \c\ 2 + \ 



2 + \e\ 2 



D k = \b\ 2 \e\ 2 + \c\ 2 \df 



"cd — bee* d* 



(10) 



Bi 



2E l 
Ck 



-A 2 
2E(A' 



A' 2 > 0, 

2 A2A 



In that case the polynomial Pk{t) (9) admits simple form 
for its roots (in particular i 2 = 0) such that we obtain 



D k = (E 2 - A 2 )(E 2 - A 



> 2 \ 



(8) 



Using relations (8), the third step consists to verify that 
the three quantities to = 4_E 2 , t\ = (A + A') 2 and i 2 = 
(A - A') 2 with < t 2 < *i < t are the three roots of 
the cubic polynomial -Pfc(t) given by 




(11) 



with A = A' and Dk > owing to i 2 = 0. The vanishing 
gap condition is thus simply 



P k (t) = * 3 - 2B k t 2 + (B 2 - 4D k )t - C 2 . (9) B 2 - AD k = (\b\ 2 + \c\ 2 - \d\ 2 - \e\ 2 ) 2 + 4\bd* + ce*\ 2 = 0. 



From this point, using standard formula (Appendix A) 
we can find the explicit form of the roots *i.2,3(fc) as 
function of coefficients Bk,C\, Dk and then obtain the 
expressions of the three quantities E, A, A' as E — ^SL, 
and A = VtT_V*i ; A' = ^E+VE f or C k > (e.g A' > A). 
The last two equalities being interchanged for Ck < 
(A' < A). We can then further substitute into (7) and 
obtain analytical but complicate expressions for energy 
bands E n (k). 

Equations (8) and (9) now constitute our starting point 
for the analysis of Dirac points properties as detailed in 
sections 3 and 4. 

In section III we consider the case of vanishing transfer 
energies along stacking j/-axis t an = (n = 1,2,3) such 
that Ck = 0. In that situation the expressions of E, A, A' 
are simple enough to derivate explicit conditions on the 
hopping transfer energies % n for the appearence, motion 
and stability of Dirac points. In the presence of inver- 
sion symmetry between A, A' molecule sites we recover 
the recent results of Mori 14 . Interestingly we obtain gen- 
eralized conditions in cases where inversion symmetry is 
absent. 

In section IV we discuss the general case of finite trans- 
fer energies along stacking j/-axis t an ^ (e.g Ck ^ 0) . In 
that situation the expressions of E, A, A' are too involved 
therefore we introduce an alternative quantity K{k) and 
derive an explicit condition on Bk , Ck , Dk for the ex- 
istence of a contact point at a position k . The main 
advantage of this condition is that it does not require ex- 
plicit diagonalization of the Bloch Hamiltonian matrix. 
Using the specific form of Bk,Ck,Dk and K(k) at the 
four time reversal points we can derive explicit conditions 
on the hopping transfer energies t an ,tb n for the appear- 
ance of Dirac point. 



III. DIRAC POINTS FOR VANISHING 
TRANSFER ENERGIES ALONG THE 
STACKING AXIS 

In this section, we consider the case of vanishing trans- 
fer energies along the stacking axis (y-axis): t an = 



(12) 

Eq. (12) admits only two kinds of solutions either (SI) 
|6| = |e| and |d| = |c| with |6| ^ \d\ or (S2) |6| = |e| = 
\d\ = \c\. A position ko, where either (SI) or (S2) is veri- 
fied, corresponds to a contact point between valence and 
conduction bands. More precisely, since the four mod- 
uli |c|, \d\, |e|) depend linearly only on k x condition 
(SI) or (S2) allows to determine the coordinate k 0x of 
position feo- The coordinate fcoj, is then determined from 



A. Dirac point condition in the presence of 
inversion symmetry between A and A' sites 

This case was recently studied in and corresponds to 
Fig. 1. The presence of an inversion symmetry between 
A and A' sites implies b = d*e lfcx and c = e *c lk:c+lky , 
(i.e., |b| = |c| = |e|) for any k. In that situation, for 
any k we obtain 



E(k) = \\b\ 



2\b\ 



2 |c| 2 - (b* 2 c 2 e- lk y + b 2 c* 2 e ik y), 

A(k) = y \b\ 2 + |c| 2 - ^2\b\ 2 \c\ 2 - {b* 2 c 2 c- lk y + 6 2 c* V fc »), 

(13) 

and the vanishing gap condition (12) simplifies to 

& = ±ieV fc */ 2 . (14) 

Equation (14) and the two equalities |6| = |c| = \e\ 
imply that in this situation the Dirac point position nec- 
essarily corresponds to a solution of type (S2) (|6| = \d = 
||c| = |e|) for (12). Furthermore, owing to the fact that 
|6|, |e| are monotonous functions of k x in [0, n] there is at 
most one pair ±fco of Dirac points such that the neces- 
sary and sufficient condition for the existence of this pair 
of Dirac points is 14 



J — JmJy < 0, 
with 

Jm = 1*6-1 ~ \tc-\ = 1*62 — *f>3 I — 1*61 — *64|, 
Jy = |*6+| — |*c+| = 1*62 + *63l ~ 1*61 + tbi\- 



(15) 
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By varying the transfer energies t bn the Dirac points 
move in the Brillouin zone as long as J < 0. From Eq. 
(14) we deduce that no Dirac point can reach X or T time 
reversal points and accordingly we also deduce that the 
case J = corresponds to a merging of the Dirac pair at 
time reversal point Y for Jy = or M for Jm — 12 ' 16 . 
The conditions Jy = or J M = define two possible 
merging pressures for each time reversal points: 



P£ = -(l±a:_)A-, 
P$ = -{l T x + )Pb+, 



(16) 



with x± given by Eq. (2) and such that J M < for < 
P < P^ { and J Y < for Py < P < P+. Depending on 
the ordering of the four merging pressures P M , P Y we 
can then deduce the pressure ranges over which the Dirac 
points existence condition J = JmJy < is verified. 
With the transfer energy values of table II, from Eq (16) 
quantitatively we obtain 



M 



-45 kbars P v 



P+ = 396 kbars 



P + 



-102 kbars, 
-15.3 kbars. 



(17) 



such that P y < P M < Py~ < P M . With this ordering of 
the four merging pressures we deduce that J = JmJy < 



for either Py < P < P M 



or Py" < P < P M . In these 



two intervals of pressure there is a Dirac pair at ±fco 
that moves by increasing pressure from Yto M for both 
intervals. Outside these intervals J > and there is a 
finite gap A between valence and conduction bands. 



B. Dirac point condition in the absence of 
inversion symmetry 

To describe a situation in which inversion symmetry is 
absent we now write 



b = t b3 + t b2 e ik * , 
c = t M c iky + t' bl e ikx+iky , 
d- 



t bi e" 
hi + h. 



°bl* 



tbl + tb4£ 



(18) 



where a priori t bn 7^ f bn for each n = 1,2,3,4. In that 
situation we need to answer if it is still possible to obtain 
Dirac points and what is or what are the necessary and 
sufficient condition(s) that generalize (15) ? 

As a first partial answer, consider the case t' bn = 
(1 + S)t bn with S ^ a real number independent of 
n = 1,2,3,4. In this situation d = (1 + 5)b*e lk * and 
c= (l + <f)e*e ifc » +ifc » (\d\ = (1 + S)\b\ and |c| = (l + 5)|e|) 
such that there is no inversion symmetry. It is immediate 
to verify that the vanishing gap condition (12) still sim- 
plifies to (14) however now the Dirac point position nec- 
essarily corresponds to a solution of type (SI) = |e| 
and \d\ — \c\ with \b\ ^ \d\). In term of transfer ener- 
gies this leads to the same existence condition (15) inde- 
pendently on 5. Thus as a first partial conclusion, this 



case shows that the absence of inversion symmetry is not 
detrimental to find Dirac point and moreover it does not 
necessarily imply more stringent existence condition. 

We now consider the general case where t' bn = (1 + 
$n)tbn with S n 7^ dependent on n = 1,2,3,4. Since 
the square moduli |6| 2 , |c| 2 , |c?| 2 , |e| 2 are linear functions 
of x — cosk x (x € [1 : —1]). to find a solution (12) of 
type (SI) (\b\ = \e\ and |d| = |c| with \b\ ^ \d\) it now 
necessitates the three following conditions: 



J -- 
J' 



(14+1 

= (K + \ 



M)(l*6-|- 

K + \W b -\ 



t z 4- t 2 
l b+ + z b 



|to-|) < 0, 

K-l) < o, 



t 2 



t 2 

l b- 



tl + + tl_ 



t', 



b+ 



li- 



ft z \fl 



(19) 

The first relation of Eq.(19) is necessary to obtain a point 
x = cosfcorr where \b\ 2 line crosses |e| 2 line = |e| at 
xq). Similarly the second relation is necessary to obtain a 
point x' = cosk' 0x where \d\ 2 line crosses |c| 2 line (\d\ = 
|c| at x' ). The last line is the necessary condition to 
obtain a unique point xq = x' . For S n (n = 1,2,3,4) 
independent on pressure, only the peculiar case S n = S 
allows to verify Eq. (19) on a finite interval of pressure 
(e.g. in that case t' bn = (1 + S)t bn such that J' = (1 + S) J 
and To = x' is automatically fulfilled). 

With the same line of reasoning, we can show that for 
generic S n dependent on n = 1,2,3,4 and independent 
on pressure, it is never possible to obtain a solution of 
type (S2) that verifies the zero gap condition Eq. (12) 
on a finite interval of pressure. 

A possible explanation for the stability of Dirac points 
against the absence of inversion symmetry is the follow- 
ing. Since transfer energies t an are ignored, the effective 
lattice model is bipartite with two sublattices composed 
on the one side by A, A' molecules and on the other side 
by B, C molecules. As a consequence of this bipartite 
property, for a nearest neighbor tight binding model there 
is necessarily a chiral symmetry S that anticommutes 
with the Hamiltonian. In our case S is represented by 
the 4x4 matrix 



S = 



( 1 \ 
10 
0-10 

\ -1 / 



(20) 



such that SH{k) + H(k)S = for t an = (e.g. a = 
f = 0). It is well established that such chiral symmetry 
favors stable band contacts. A confirmation of the role 
played by this chiral symmetry would be to examine the 
existence and stability of Dirac points when a small onsite 
potential is added with (Va = Va') (Vb = Vc) since 
such a perturbation commutes with chiral symmetry and 
furthermore annihilates the bipartite property. 
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IV. DIRAC POINTS FOR FINITE TRANSFER 
ENERGIES ALONG THE STACKING AXIS 



for p — 0, 1, 2 and where 



We now examine the general case of finite transfer en- 
ergies along stacking y-axis t an ^ (n = 1,2,3). (e.g 
a^O and / ^ such that Ck ^ 0). In that situation, 
using so called Viete-Descartes formula (see Appendix A) 
we can obtain the roots of cubic polynomial Eq. (9): 



t P (k) = 3^ + 



Bl 



12Dk cos ( 



k - 2pn 



)], (21) 



cos Ok 



TlCj + 11D k Bk - 2Bj 
2{B 2 k + 12^)3/2 ' 



(22) 



For Cfc > (A' > A) we can then deduce explicit ana- 
lytical formula for E, A, A' valid for any k: 



E(k) 
A(k) 
A'(fe) 



Via _ / Bfc + \/Bf+l2gIcos(%) 
2 V 6 

Jtj-Jh. _ / j3 fc +^B£+i2n fc co S (' 

2 V , 6 



VtT+V^ _ B h + ^/Bl + 12D h co S (-K 
2 V 6 



B k +y/Bl + 12D h cos(g*4 

6 

Bfc+V-Bg + 12-Dfc cos(^ 



(23) 



The last two equalities being interchanged for Ck < 
(A' < A). We emphasize that Eq. (23) can be further 
substituted into Eq. (7) to obtain analytical expressions 
for energy bands E n (k). 



A. Contact point condition 

Interestingly, from Eq. (23) we obtain that the condi- 
tion of vanishing gap A = {E\ = Eq) is Ok = and 
Ck > 0. Similarly the condition A' = (E 3 = E4) is 
Ok = and Ck < and the condition E2 = E3 is Ok = 7r. 
Starting from Eq. (22) it is easily shown that we can 
rewrite (B 2 k + 12L> fc ) 3/2 (l - cos6>) = K+(k)K-(k) such 
that K±(k) > vanishes when Ok = and TCfc > 0. 
For our main focus being contact points between va- 
lence and conduction bands we only consider K_(k), 
such that conditions 0k — (Ck > 0) can be written 
as K(k) ee K (k) = with 




FIG. 2. Parametric curve yk{xk) solution of K(k) — 0. From 
Eq. (8) in the case A = we obtain the parametric forms 

X W = 4 g (1 + t/ t 2 )3/2 ,yW = JTWW 1 G P : 4 

t = A' 2 /E 2 . 



K{k) = B 3 k /2 [(2 - A /TT3^)Vl + VT+Wk 



where yk 



4D k 



(-§ < Vk < 1) and x k 



2x k ] , 
_(24) 

27 C fc 

8 R 3/2 



(0 < a;^ < 1) (in Appendix B we present an alterna- 
tive derivation of the quantity K(k) 19 ). It is immedi- 
ate to verify that for Ck = (a;*, = 0) the condition 
K(k) — is equivalent to B k = 4Dk {yk = 1) which 
is exactly Eq. (12) found in section 3. Other simple 
cases that verify K(k) = are Xk = l,yk = —1/3 and 
Xk = I/a/2, Vk = 0; as we show below this last case is re- 
alized at r and X for t a \ =0. More generally the values 
Xk,Vk such that K(k) = define the parametric curve 

x(t) = (1+t / 2) 3/2 , y(t) = j^tJW for * e t° : 4 1 shown 
on Fig. 2. The main qualitative feature to retain from 
this curve is that yk is a monotonous decreasing function 
of Xk {Ck) moreover an approximate linear interpolation 



4-Dfc 
Bl 



1 



m c k 



B 



3/2 , 



More 



gives y fe ~ 1 - V2x k (e.g. 

quantitatively, as we show in next sections, this para- 
metric constraint allows to derive explicit conditions in 
terms of t an , t on for the merging conditions at M and Y 
points and similarly at T and X for t a \ = 0. In practice, 
as illustrated by Fig. 3 obtained from transfer energies 
tan, tbn of table I at P = 4 kbar, we can locate Dirac 
points position ±fc by locating the positions of zeros of 
K(k) > without direct diagonalization of the Bloch 
Hamiltonian 4x4 matrix. We note that in the neigh- 
borhood of a Dirac points, we have 0k <C 1 such that 
K(k) oc (1 — cos^fc) oc 0\ and A(fc) cx smOk oc 0k (from 
Eq. (23)), thus A(fe) cx y/K{k). 
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FIG. 3. (Color online) Contour plot of J K(k) obtained from 
transfer energies t an ,t^ n of table I at P = 4 kbar and for 



y in the first Brillouin zone (— n < 



k x , k 

center corresponds to the F point. 



k x , k y < 7r). Fhe 



B. Dirac point in the presence of inversion 
symmetry between A and A' sites 

1. Merging condition and merging pressure at Y and M 
points 

We have seen in section III A that Dirac points only ap- 
pears over finite intervals of pressure. In the presence of 
time reversal and inversion symmetries the lowest (high- 
est) pressure of each interval necessarily corresponds to 
the merging (or emerging) of a Dirac pair at a time rever- 
sal G/2 point of the Brillouin zone 16 . It is thus essential 
to determine the conditions for merging at each of the 
four time reversal points G/2 — V, X,Y, M; since the 
possibility of having Dirac points in the Brillouin zone 
depends on the existence of such merging points. 

For the case t an = we have seen in section III A that 
only Y and M time reversal points constitute merging 
points. Moreover at such merging points the Dirac point 
condition Eq. (12) becomes the merging conditions Jy — 
|*6+| - |*c+| = at Y and J M = |i 6 _| - \t c -\ = at M 
Eq. (15). We now explain how these merging conditions 
are modified for t an ^ 0. Since the reasoning is similar 
for M and Y we describe it only for M point. To start 
with we note that at M we have / = 0, b = —d = 
c = e = < c _ and a — t a —. In that situation the three 
quantities E>m,Cm,Dm read: 



B M =a 2 + 2(b 2 +e 2 ), 
C M = 2a(b 2 -e 2 ), 
D A4 = Ae 2 b 2 , 



(25) 



from which we obtain 



Vm 
xm 



AD M _ 4e 2 b 2 
B M (^#+6 2 +e 2 ) 2 ' 
V27 C M _ %/27 a(b 2 



B 



4 (^+fc 2 +e 2 ) 3 / 2 ' 



From these expressions by noting t = 



fc 2 +e 2 
(h 2 -e 2 ) 



vided that a verifies the equality a 
recover the parametric form xm = —4— ( 1+t /2) 3 / 2 



(26) 



and pro- 
then we 



Vb 2 +e 2 

^ ' and 



l-t 



Vm — (i+t/2) 2 



which is equivalent to Km = 0. We thus 



(b —e ) 

conclude that the merging condition at M is a = \j b2+ J 2 ■ 
For Y point, with a similar reasoning we find a merg- 
ing condition a = — ^p=f with b = t^+, e = t c+ and 
a = t a ~. 

In summary, for t an ^ we obtain modified merging 
conditions Jm — 0, Jy = with: 



Jm — t 2 _ — t 2 _ — t a - 
Jy = t b+ - t 2 c+ + t a _ 



(27) 



For < a 3 = t a 2 = we recover the condition Eq. (15). 
From Eq. (27) it is possible to obtain approximate ana- 
lytical expressions of the merging pressures Pm, Py that 
verify J M {P M ) = and J y {Py) = 0. The detailed 
derivation is described in the Appendix C. At M point, 



we distinguish two cases: (Ml) for 



< 1 and 



(M2) for 



C- n- 



» 1. The latter case (M2) appears 



in line with parameters given in table II, moreover since 



Pa- 



Clwc obtain 



l-xt 



P, 



M 



\P h -\ 



(28) 



At Y point, we also distinguish two cases: (Yl) for 



C- P b + 
t° + P a - 



< 1 and (Y2) for 



P b + 
P a - 



1. Again the 



latter case (Y2) appears in line with parameters given in 



table II and moreover since 



Pa- 



C 1 we obtain 



\P h - 



v 71 



(29) 



Expressions (28,29) provide a non trivial dependency of 
the merging pressures Py,Pm in terms of the parame- 
ters (tP a _ 1 P a _) that characterize tranfer energies along 
the stacking axis. Introducing a rescaling parameter r 



such that t° a _ 
to P M Py 



rta_, for large \r\ both expressions lead 
-P a -. For r = + (such that the denomi- 



P, 



M- 



nators never vanish) Eqs.(28,29) lead to Pm — -^r 
and P Y = ^ ± P Y+ - Fi gures 4a and 4b further illustrate 




read 
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FIG. 4. (Color online) r is a rescaling parameter r such that 
t°_ — > rt° a _. (a) grey and white regions indicate respectively 
the regions Jm(P,t) < and Jm(P,t) > obtained from 
Eq.(27). The curves P M (r) (blue line), obtained from Eq.(28), 
accurately define the frontiers separating Jm < and Jm > 
regions, (b) Similar diagram at Y point. Regions Jy < and 
Jy > are obtained from Eq.(27) and the curve Py{t) (red 
dashed line) are given by Eq.(29). 



the validity of Eq. (28,29). On Fig. 4a, in a diagram 
(P, r), the regions Jm(P,t) < and Jm{P, r) > ob- 
tained from Eq.(27) are indicated in grey and white re- 
spectively. On the same Fig. 4a the curves Pm{t) (blue 
lines), obtained from Eq.(28), accurately define the fron- 
tiers separating Jm < and Jm > regions. Fig. 4b 
show similar results at Y point. 



2. Merging condition and merging pressure at V and X 
points 

At r and X points we have / = 2t a i and a = t a+ . As in 
the previous case we consider only T since the reasoning is 
similar for X point. At T, the three quantities -Br, Cr, P>t 



B r = f 2 + a 2 + 2(b 2 + e 2 ), 
C r = -2(2feb + a(e 2 + b 2 )), 
Dv = fa(fa-Aeb), 



(30) 



with e = t r 



Vv 



tb + . We then obtain 

_ /a(/q-4eb) 



y/27 2/efc+q(e 2 +b 2 ) 
4 (^^+b 2 +e 2 ) 3 / 2 ' 



f 2 +a 2 



(31) 



then 



b 2 +e 

such that 



In the spirit of previous section we define t 
the aim is to find a condition Jr(f, a, b, e) 
we can recover the parametric form yr — y{t)i %t = 
x(t). For t a i ^ 0, we could not find such self consistent 
solution explicitly. However, for t a i = (/ = 0) we 
immediately obtain that the self consistent solution is 
t = 1 (e.g. dy = 0, xr = I/n/2) which necessitates the 
condition Jp = \Jb 2 + e 2 + a = 0. In a similar manner 
at X point we find the condition Jx — yjb 2 + e 2 — a = 
with e = — t c _, b = t(,_. In summary for i al = 0, we 
obtain the merging conditions Jr = 0, Jx = with: 



for t a i = 



J) 

Jx 



r - \l t l+ + f 'c^ 



t„ 



4- + tl- 



(32) 



a+ 



From Eq. (32) it is possible to obtain approximate ana- 
lytical expressions of the merging pressures Pp , Px that 
verify Jr(-P r ) = °> Jx(Px) = 0. At T (X) we find that 



the condition 



>1( 



> 1) is necessary. 



These two inequalities are indeed verified by parameters 
given in table II. Assuming further that \Pr/Pb+\ *C 1 
and \Px/Pb- \ <C 1 we then obtain : 



for t a i =0 



Pt 



1 



1*2+1 



P 



X 



+ 



'1*2+1 Pa+ 
t° 

_ °+ 

I'M 
*2+ jn^i 

i*o p a + 



P, 



P h -\ 



(33) 



Equalities (33) provide a non trivial dependency of the 
merging pressures Pr — Px in terms of the parameters 
(t a+ ,P a +) for the case t a \ = 0. Introducing a rescal- 
ing parameter r such that — > rt° a+ , for large \r\ 
both expressions lead to Pr — Px — —Pa+- Figures 
5a and 5b further illustrate the validity of Eq. (33). On 
Fig. 5a, in a diagram (P, r), the regions Jr(P,r) < 
and Jr(P, r) > obtained from Eq.(32) are indicated 
in grey and white respectively. On the same Fig. 5a 
the curves Pr(r) (magenta dotted line), obtained from 
Eq.(33), accurately define the frontiers separating Jr < 
and Jr > regions. When t a \ ^ 0, Eqs. (32,33) are no 
more valid and the merging pressure Pv(r) is directly 
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of transfer energies t a ± the following four inequalities 



FIG. 5. (Color online) r is a rescaling parameter r such that 



t"_i_ — > rt^_|_. (a) when t a i = : grey and white regions in- 



dicate respectively the regions Jr(P, r) < and Jr(P,r) > 
obtained from Eq.(32). The curves Pr(r) (magenta dotted 
line), obtained from Eq.(33), accurately define the frontiers 
separating Jr < and Jr > regions. When t a i 7^ : ma- 
genta x and magenta + show the merging pressures Pr(r) 
calculated numerically from the condition Kr (P, r) — for 
the cases ±t a i with t a i taken from table I. (b) Similar dia- 
gram at X point. 



calculated from the condition Kr(P,r) = 0. The corre- 
sponding numerical results Pr(r) for ±i i (magenta x 
and magenta +), with t a \ taken from table I, are shown 
on Fig. 5a for comparison. Fig. 5b show similar results 
at X point. These figures show that by changing the sign 
of r, the role played by T and X points are interchanged; 
more precisely Px < —P a + < Pr when r > 0, while 
Pr < -P a + < Px when r < 0. 



Before going further we summarize the main finding of 
the last two sections. We have seen that in the presence 



at M: 
at Y: 
at T: 
at X: 



,t *- n- 

I t°_ Pa- 
1*2- P b+ 



1 t° b+ Pa+ 

\wzp^ 



> l, 

> l, 

> l, 

> l 



(34) 



constitute the necessary conditions to find at least and 
at most one solution to the merging conditions Jg/2 = 
at the four time reversal points M, Y, T and X. Beside 
these necessary conditions, the sign of t° a+ plays an im- 
portant role. For t° a+ < as given in table II, it appears 
that only Pr can take a positive value (physically acces- 
sible pressure). For t° a+ > 0, there is a small region of 
parameter where both Px,y,m can take positive values. 
Finally we note that transfer energy t a i seems to play 
only a marginal role in the determination of the merging 
pressures Pp,Px (t a i plays no role in the determination 
of P y ,m)- 



3. Existence condition and existence domain of Dirac 
points: 

In the absence of transfer energies t an we have seen 
that the condition Eq. (15) for the existence of Dirac 
point in the Brillouin zone is obtained as J = JmJy < 0. 
In the presence of non zero transfer energies t an a direct 
generalization of this condition reads 20 



J — JmJyJxJt < 0, 



(35) 



with Jq/2 the quantities defined by Eq.(27) and Eq.(32) 
in the preceding sections. For t an = 0, Eq. (32) show 
that Jr.x > such that the condition Eq.(35) becomes 
equivalent to the condition Eq. (15). 

More generally, for t a \ — 0, in the spirit of the pre- 
vious sections Fig. 6 shows the regions J(P, r) < (in 
grey) and the regions J(P, r) > (in white) obtained 
from using J = JmJyJxJt] with r the rescaling param- 
eter defined before (t° ± — > rt° a± ). On the same Fig. 6, 
the merging pressures curves Pg/2{t) deduced from Eqs. 
(28,29,33) are also plotted. As espected, the merging 
pressures curves define the frontiers separating regions 
J < from regions J > 0. Again, we emphasize that in 
regions J < there is a pair of Dirac points in the Bril- 
louin zone whereas in regions J > there is a (direct) 
gap separating valence and conduction bands at any k. 
More concretely, as an example of interpretation of Fig. 
6, when r = 1 the region J < corresponds to the pres- 
sure intervals Px < P < Pm and Py < P < Pr', as a 
consequence there is pair of Dirac points (in the Brillouin 
zone) in these two pressures intervals. More precisely, 
upon increasing pressure, there is a pair of Dirac points 
emerging at X (Y) at Px (Py) and merging at M (T) at 
Pm {Pr)- When r > 2, Fig. 6 shows that the ordering of 
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FIG. 6. (Color online) t a i = 0: diagram (P, r) showing the 
regions J(P, r) < (in grey) and the regions J(P, r) > 
(in white) obtained from using J = JmJyJxJt', with r the 
rescaling parameter defined before (t°± — > fta±)- The merg- 
ing pressures curves Pg/2(0 deduced from Eqs. (28,29,33) 
are also plotted. As espected, the merging pressures curves 
define the frontiers separating regions J < from regions 
J > 0. 



the merging pressures has changed such that the region 
J < corresponds now to Py ps P m and Px < P < Pr- 
This implies that the motion of Dirac points in the Bril- 
louin zone is strongly affected by the explicit value of r 
(e.g. t°±) all other parameters being kept fixed (see also 
the case of negative r). From a more physical point of 
view, Fig. 6 shows that the condition r > is more favor- 
able to observe Dirac points, since the region J < has 
a greater overlap with the region of physically accessible 
(positive) pressure values. 



C. Stability of Dirac points in absence of inversion 
symmetry 

In the absence of transfer energies t an we have seen 
that the loss of inversion symmetry was not detrimen- 
tal to the existence of Dirac points and we attribute this 
stability to the presence of a chiral symmetry. In the 
presence of finite t a ± this chiral symmetry is lost because 
the corresponding nearest neighbor tight-binding model 
is no more bipartite. As a consequence the absence of 
inversion symmetry might now be sufficient to prevent 
the existence of Dirac points. We note however that the 
possibility of Dirac points in absence of both inversion 
and chiral symmetries was recently exemplified for a two 
bands tight binding model on the Honeycomb lattice 21 ; 
according to the authors the stability of Dirac points re- 
lies on the existence of an average-inversion symmetry. 
We now explore this possibility in our four bands model. 

As in section III B to describe a situation in which 
inversion symmetry is absent we now rewrite (b, c, d, e) 
as in Eq. (18) such that \b\ ^ \d\ and |c| ^ |e| for 



all k. Furthermore we also define b + = b(k x = 0) and 
6_ = b(k x = 7r) and similarly we define c±,d±,e±. We 
now reexamine the merging properties at M. The three 
quantities Bm,Cm,Dm read 



B M = a 2 _ + b 2 _ + c 2 _ + d 2 _ + e 2 _, 
Cm = 2a_(e_c_ — b-d-), 
D M = (e_b_ +c_oL) 2 , 



(36) 



from which we obtain 

4(e_b_+c_d_) 



yM — - ( a 2_ +i ,2_ +c 2_ +d 2 +f ,2J2 , 

Tf C M _ v / 27 o3/2 a_(e_c_-6_d_) 
8 — 4 Z (a 2 _ +b 2 _ +c 2 _+d 2 _ +e 2 _ ) 3 / 2 ' 



(37) 



M 



From this point by noting t 



2a 2 _ 

b 2 +c 2 ~+d 2 +e 2 



if a_ 



y / 2(e_c_-b_d_) 

^ b 2_ +c 2 +d 2_ +e 2_ 

Vr? t 



and if b 2 _ + c 2 _ — d 2 _ + e 2 . then we 



can rewrite x M = ^j 1 n+hm and Vm = (i+W which 
implies Km = 0. We can proceed similarly for the other 
time reversal points such that when inversion symmetry 
is lost the generalized merging conditions read (for T and 
X point we take t a \ — 0) 



at M: Jm — (e_c_ — b-d-) — fl- 
at Y: Jy = (e+c+ — b + d + ) — a^^Jb\ + c 2 
at V: Jr = (e + c + + b + d + ) + a + Jb 2 + + c\, 
at X: Jx — (e_c_ + b-d-) + 



(38) 



where we have taken care of the two supplementary con- 
straints 



A- 


+ c 2 _ 


k\- 




b\ -\ 


-4. 


lb 2 _ 





bl 



(39) 



Assuming that these two constraints are verfied for any 
pressure it can be rewritten as a single constraint valid 
for any k x 



\b\' 2 



\df 



(40) 



Summarizing, we found that the loss of inversion symme- 
try does not prevent the existence of Dirac points pro- 
vided a kind of average-inversion symmetry Eq. (40) is 
verified; this is reminiscent of what was recently obtained 
on a two band model on the honeycomb lattice 21 . As a 
last caveat, we emphasize that it is possible to construct 
bands model (time reversal symmetric) but with neither 
inversion nor chiral symmetries, and that exhibit acci- 
dental but stable Dirac points (with zero gap state) that 
cannot merge at time reversal points. 



V. SUMMARY AND DISCUSSION 

In this work we developed a new method to study the 
existence, stability and merging of Dirac points between 
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valence and conduction bands of 3/4 filled a-(BEDT- 
TTF)2l3 conducting plane. We considered the usual 
nearest neighbor tight binding model with the seven 
transfer energies t an (n = 1, 2, 3) ,tb n (n = 1, 2, 3, 4) that 
depend on the applied pressure with associated charac- 
teristic pressures P an ,Pbn- Owing to the four distinct 
molecules (A, A', B, C) per unit cell of the Bravais lat- 
tice, the corresponding Bloch Hamiltonian is a 4 x 4 ma- 
trix H(k) for each wave vector k of the Brillouin zone. 
In most previous works the study of the Dirac points was 
achieved through direct numerical diagonalization of this 
matrix "H(fc). In this work we have shown that analyti- 
cal understanding of the physics of Dirac points is within 
our grasp. In a first part we have thus explained how it 
is possible to obtain formally complete analytical albeit 
complicate expressions of the four energy bands. 

As a first application of our method we reexamined 
the simple case where transfer energies along the stack- 
ing axis are ignored t an = 0. In the presence of inversion 
symmetry we recovered the Dirac points existence con- 
dition Eq. (15) recently obtained by Mori 14 . We have 
shown that this condition Eq. (15) defines the existence 
domain of Dirac points from the merging conditions at 
M and Y time reversal points. We then considered sit- 
uations in which inversion symmetry is lost due to an 
increased degree of anisotropy between the transfer en- 
ergies tbn- We have shown that the absence of inversion 
symmetry is not detrimental to the existence of Dirac 
points and we derived a generalized existence and stabil- 
ity conditions Eq. (19). A possible explanation of the 
stability of Dirac points in the absence of inversion is the 
existence of a chiral symmetry, the latter being present 
owing to the bipartite property of the system when trans- 
fer energies t an vanish. This idea needs however to be 
further explored. 

In a second and main part of this work we considered 
the general situation with the seven transfer energies. 
The analytical expression of the gap A(fc) between va- 
lence and conduction bands being too involved to analyze 
the existence of Dirac points, we proposed an alternative 
quantity K{k) > Eq. (24) that vanishes only at the po- 
sition of band touching points between valenc e and con- 
duction bands. More quantitatively A(fe) oc \J K(k) near 
band touching points. Using this alternative quantity, in 
the case where inversion symmetry is present, we deter- 
mined the merging conditions at each four time reversal 
points since knowning the merging properties at the time 
reversal points is sufficient to determine the existence do- 
main of Dirac points 20 . For M and Y points we obtained 
generalized merging conditions as compared to the first 
part Eq.(27) vs Eq.(15). From these merging conditions, 
we derived analytical formula for the associated merg- 
ing pressures Pm Eq. (28) and P Y Eq. (29) functions 
of the different transfer energy parameters t an ,tb„ and 
their associated characteristic pressures P an ,Pbn- For the 
case t a \ = 0, we obtained merging conditions at Y and 
X Eq. (32) and derived analytical formula for merging 
pressures Pr and P a Eq. (33). As exemplified by Figs. 



(4a,4b,5a,5b), the analytical formula Eq. (28,29) and Eq. 
(33) (for t a \ — 0) appeared to agree perfectly with numer- 
ical results. All in all, from the analysis of the merging 
conditions and merging pressures at each time reversal 
point we deduced four inequalities constraints Eq. (34) as 
necessary conditions. We also emphasize the importance 
of the sign of t a ± in the value taken by the merging pres- 
sures; by contrast we also pointed out the marginal role of 
transfer energy t a \ on the merging pressures. Combining 
the merging conditions at the four time reversal points, 
we proposed a generalized Dirac points existence condi- 
tion Eq.(35) as compared to Eq. (15). This condition 
Eq.(35) is reminiscent of the recent results obtained in 20 
with a totally different method. An illustration of this 
existence condition is given by Fig. 6 (for i al = 0). We 
explained how from Fig. 6 it is possible to determine the 
intervals of pressure allowing for the presence of Dirac 
points in the Brilloin zone. We then exemplified how the 
allowed intervals of pressure depend on the ordering of 
the four merging pressures. We stress that a more thor- 
ough exploration is needed to study the motion of Dirac 
points in these allowed pressure intervals. In particular it 
would be important to discriminate interval of pressure in 
which Dirac points are at the Fermi level (so called zero 
gap state) from pressure range in which Dirac points be- 
tween valence and conduction bands exist but stay below 
(above) the Fermi level (metallic phase with electron-hole 
pockets). At last we also explored the role of inversion 
symmetry in the presence of transfer energy t an . We 
have shown that the loss of inversion symmetry does not 
prevent the existence of Dirac points provided a kind of 
average-inversion symmetry is maintained Eq. (40). In 
that condition we have established generalized merging 
conditions at each time reversal point Eq. (38). A more 
thorough study is needed to understand the opening of 
the gap in case where the average invcrision symmetry is 
not fulfilled Eq. (40). More generally a natural extension 
of this work would be to consider the effect of anions 22 
induced onsite potentials as well as mean field interaction 
effect. 
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Appendix A: Viete-Descartes roots formula for 
cubic polynomial with three real roots 

We briefly reminds the derivation of Viete-Descartes 
formula for the roots of a cubic equation 



t 3 + at 2 +bt + c = 0. 
First we change the variable from t to r\ as 

t 





with p = a 2 /3 — b. The equation reads 
Arj 3 - 3r? 

with q = -2a 3 jll + ab/3 - c. For 
introduce 9 by 

COS0 

and then rewrite 

4?7 3 - 3?7 = 4cos 3 

The roots of this equation are 

9 (9 2vr\ (6 4tt 

7] = cos - , cos — , cos — 

' 3' \3 3/' \3 3 



(Al) 



(A2) 



(A3) 



^q\ < 1 we can 




- — 3 cos -. 
3 3 



(A4) 



(A5) 



(A6) 



In our case we have a = —2B kl b — B\ — 4Z? fc and c = 
-C\ such that p = and q = ^g£±^^gl , 

Eq.(A4) is then equivalent to Eq.(22). 



Appendix B: Alternative derivation of quantity K(k) 
Eq. (24) 

In this appendix wc present another derivation 19 of 
Eq.(24) that does not necessitate the use of cubic polyno- 
mial Eq.(9) . We start from the characteristic polynomial 
of the 4x4 Bloch hamiltonian matrix: 

Fk(uj) = lu a - B k uj 2 + C k Lo + D k 

= (w - E 1 {k)){u - E 2 (k)){uj - E 3 (k))(uj - E 4 (m) 

We reminds that B k > and we restrict to the case 
Cfc > which is a necessary condition to obtain a contact 
point between valence band E 2 and conduction band E\ . 
We now define three functions: 

Fi M = = 4w 3 _ 2BkUJ + Ck ; (B2 ) 

*fcM = F k(^ - F h (w) = 3uj 4 - B k u? - D k , (B3) 
Fh(u) = F k(")u - 4 ^fcM = 2Bfc^ 2 - 3C k uj - 4^4) 



We denote loq, u a and uj], the largest root of each of 
these functions respectively (Fi(u>o) — 0, F a (uj a ) = 
and F b (u) b ) = 0). Quantitatively we obtain 



Bk + y/B\ + 12£> fc 

6 ' 

3C fc + y/9C'j + 32B k D k 
45 fc 



(B5) 
(B6) 



By construction we have E 2 < ujq < E\ since F k (E 2 ) < 
and Fl(E x ) > 0. From F k (E 2 ) = and F k (uj ) < 
we further deduce that F^(E 2 ) < and F£(w ) > 
implying E 2 < uj a < luq < E±. Very similarly we further 
obtain E 2 < uj b < cu a < lu < E\ owing to F k (E 2 ) < 
and F^(u! a ) > . Wc thus deduce that if E x — E 2 then 
necessarily u>b = ui a - By defining 19 



k(k) =v /, 



1 + VI + 3?7fc - {xk + V x l + 3yfe) 2 

(B7) 

with y k = and x k = y^^pyj ■ We then obtain that 

K(k) = when E\ = E 2 . It is straightforward to show 
that K(k) = is equivalent to K(k) = Eq.(24). This 
derivation appears simpler than the one presented in the 
bulk of the article however it remains a caveat which is 
to demonstrate that E\ — E 2 when K(k) = 0. 



Appendix C: Explicit determination of merging 
pressure at M point 

In this appendix we present the detailed derivation of 
merging pressure Pm (Eq.(28)) that verifies Jm(Pm) = 
with 



J M {P)=tl_-tl_-t a .^tl_+tl 



(CI) 



(see Eq. (27) of the main text), where t a _(P) = 
*2_ (1 + pB' = + Tt) and t c _(P) = 

t° c _ with t° al P a a = (a_,&_,c_) given in table II. 
Defining x = P/\Pb-\, we rewrite Jm(P) as Jm{x) = 
(tDHgix) - f{x)) with g{x) = (1 - x) 2 - x 2 _ and 

/(*) = + ^x)^J{l-xf +x 2 _. A pressure P M 

solution of the merging condition Jm(Pm) — then cor- 
responds to a point xm = PAi/\Pb-\ where the curve 
g(x) intersects the curve f(x). To determine the posi- 
tion xm we note the following properties of functions 
g(x) and f(x). The function g(x) is a parabola that 
has two zeroes at xm± = PM±/\Pb-\ = — (1 ± x -), a 
minimum at x = 1 and g(x) cx x 2 for large |x|. The 
function f(x) is monotonous, it has a single zero at 

x a - = -r§f| a nd f(x) cx ^T\^T X \ X \ for lar § e l x l- 
With the parameters as given in table II the zeroes xm± 
and .x a _ verify — 1 < Xm- < x a _ < 1 < x m+- From 
these properties we can formally distinguish two cases: 
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(Ml) for 



\Pb-\ 



< 1 and (M2) for 



ta- in- 1 



» 1, 



|*& — I Pa- ^ V ; 1*6-1 Pa 

only this latter case appears in line with a-(BEDT- 
TTF) 2 l3 parameters as given in table II. In the relevant 
case (M2) the condition Jm(%) = has only one solu- 
tion xm- For i a _P a _ > (M2a) this solution verifies 
— 1 < Xm- < x M < x a _ < 1 whereas for t a _-P a _ < 
(M2b) one finds x a - < Xm < %m+- Strictly speaking 
only (M2a) corresponds to a-(BEDT-TTF) 2 I 3 parame- 
ters as given in table II; in that situation since \xm\ < 1 
we can linearize J M (x) to obtain an approximate expres- 



sion of the merging pressure Pm- 



M 



2+7 



\Pb-\ 'a- 



(C2) 



This expression corresponds to Eq.(28) given in the main 
text. For the Y, T and X points we can proceed very 
similarly from their respective condition Jy,v,x{P) = 
and obtain the merging pressures formula Pm , Pr and 
Px Eq.(29) and Eq.(33) given in the main text. 
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